Mottness scenario for non-Fermi liquid behavior in the Periodic Anderson Model 

within Dynamical Mean Field Theory 
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We study the Mott metal-insulator transition in the Periodic Anderson Model within Dynamical 
Mean Field Theory (DMFT). Near the quantum transition, we find a non- Fermi liquid metallic state 
down to a vanishing temperature scale. We identify the origin of the non-Fermi liquid behavior as 
due to magnetic scattering of the doped carriers by the localized moments. The non-Fermi liquid 
state can be tuned by either doping or external magnetic field. Our results show that the coupling 
to spatial magnetic fluctuations (absent in DMFT) is not a prerequisite to realize a non-Fermi liquid 
scenario for heavy fermion systems. 
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The theoretical understanding of the breakdown of the 
Fermi liquid paradigm observed in high Tc superconduc- 
tors and heavy fermions systems remains one of the open 
challenges in strongly correlated physics. These systems 
show metallic phases with anomalous properties that can- 
not be accounted for by the Fermi liquid theory, which 
provides an adequate description of the electronic state 
of ordinary metals. The reason is intimately related to 
the strong correlation effects, originated in the localized 
nature of the d and / orbitals of the experimental com- 
pounds [ij. Different ideas have been proposed over the 
years to try to explain the origin of the non-Fermi liquid 
(NFL) states, without an absolute consensus so far. How- 
ever, many of these ideas share a common feature, namely 
that the central ingredient is the proximity to a quan- 
tum phase transition (QPT), or a quantum critical point 
(QCP) ^2]. In that scenario the breakdown of the Fermi 
liquid occurs in the neighborhood of a T = transition 
between an ordered phase {e.g. antiferromagnetic) and 
a paramagnetic one. There, the fluctuations of the or- 
der parameter that couple to the electrons are strongest, 
and are viewed as the origin for the NFL state. Among 
those approaches we can mention the Hertz-Millis the- 
ory, where the paramagnons of the ordered phase "dress" 
the conduction electrons to produce the NFL features 
0. Another approach is the local quantum critical the- 
ory 0, [sj , which does not consider the electrons as "by- 
standers" but emphasizes their role in the screening of 
the local magnetic moments, via the celebrated Kondo 
effect. There, the competition between the tendency to 
formation of local singlets and the long wavelength mag- 
netic fluctuations are to be considered on equal footing. 
This has been achieved by the formulation of an exten- 
sion to the Dynamical Mean Field Theory (DMFT) ap- 
proach called EDMFT 0, H- DMFT has proven to 
be a very useful technique to study strongly correlated 
electron systems when the main physical effects are local 



jH]. However, it is also recognized that this method lacks 
a proper description of the spatial magnetic fluctuations 
which are usually considered a crucial ingredient for the 
realization of a NFL state. This shortcoming is cured 
in EDMFT with the incorporation of a bosonic compo- 
nent to the effective mean field. Though that approach 
has provided useful insight into the physics of the prob- 
lem 0, H, 1^ , the solution of its mean field equations 
is more difficult and often requires some simplifying as- 
sumptions. We should also mention that a different and 
original approach to quantum criticality that goes beyond 
the Ginzburg-Landau theory has also been proposed [lH . 

In the present work we shall show that, contrary to 
conventional expectations, a non- Fermi liquid state is 
readily obtained from the DMFT solution of the canoni- 
cal model for the study of heavy fermion systems, namely, 
the periodic Anderson model (PAM) [12|] . We shall show 
that, similarly to other approaches this novel NFL state 
is located in the neighborhood of a QPT, but unlike the 
standard scenario described before, the relevant quantum 
transition here is a Mott transition. Thus, the present 
study sheds a different light onto the problem, showing 
that the coupling to long wavelength magnetic fluctu- 
ations (which are absent in DMFT) is not a prerequi- 
site for the realization of a NFL scenario that captures 
some features observed in heavy fermion systems. Local 
temporal magnetic fluctuations alone can provide sufB- 
cient scattering to produce an incoherent metallic state. 
Two recent works have also reported non-Fermi liquid 
states in models that generalize the PAM by introduc- 
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ing a finite bandwidth to the localized orbitals 
However, those approaches fundamentally differ from the 
present one, since they view the NFL state as the result 
of an orbital selective Mott transition, which is absent in 
the PAM. Thus, our study shows that the PAM, solved 
within DMFT {i.e. in the limit of infinite dimensions), 
may be considered as a "bare bones" or minimal ap- 



proach that can qualitatively capture a physical scenario 
found in some NFL heavy fermion systems. 

The periodic Anderson model Hamiltonian reads, 



H = -t {ptpja + PjaPia) + (cp - fi)J2ptP^'y 

<ij>(7 iCT 

ia 1(7 

i 

(1) 

where and are the destruction and creation oper- 
ators for the electrons at p orbitals with energy Cp and 
hopping parameter t. dia and are the respective oper- 
ators of the non-dispersive c?-electron orbitals with en ergy 
Cd, that we fix equal to without loss of generality [15l |. 
The p and d orbitals are hybridized with an amplitude 
tpdi and the electron correlations are introduced by the 
Coulomb interaction U on the d sites only. For simplic- 
ity we consider the model defined on a Bethe lattice that 
gives a semicircular DOS for the p-electron band (at tpd 
and U = 0) [6i] . Its half-bandwidth is D = 2t = \ and sets 
the units of the problem. The charge transfer energy is 
defined as A = — fp, and /x is the chemical potential. 
At [/ = and large A, the hybridization parameter tpd 
permits the delocalization of the d electrons, as a narrow 
band forms at the Fermi energy with bandwidth ~ tp^/A. 
However, when both U and A are large, and the narrow 
band is half-filled, the periodic Anderson model can de- 
scribe two different correlated insulator states: a charge 
transfer insulator for A < [/, and a Mott-Hubbard one 
for U < A. 

The doping driven metal-insulator transition (MIT) in 
the paramagnetic Mott-Hubbard insulator was the focus 
of our recent study . The main finding was a quali- 
tatively different scenario for the electron or hole driven 
transitions. In the former case the MIT was expectedly 
similar to the first order transition of the well studied one 
band Hubbard model Q. However, in the latter case, an 
intriguing second order transition was found. Here, we 
shall study in detail the unexpected nature of this Mott 
transition. 

Within DMFT the PAM can be mapped onto a quan- 
tum inipurity problem subject to a self-consistency con- 
dition [l7|. Using the cavity method Q one obtains the 
effective action for the quantum impurity model: 



r0 r0 

Seff = - / dr dr' J2 dt {r)G^\r - T')d, (r') 
Jo Jo ^ 

+ U f dT K(T)-i] [ndi{T)-\] 
Jo 



(2) 



where d^ and destroy and create a d electron on an 
arbitrarily chosen site and Uda- is the occupation opera- 
tor. Qq is the function describing the properties of the 
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FIG. 1: Density of states for p and d electrons (solid and 
dashed line respectively) from QMC data at T=l/128. The 
doping is 5 = 0.049, fi = 0.53, A = 1, [/ = 2 and tpd = 0.9. 
The analytic continuation is done using maximum entropy 
method. Inset: DOS of p and d electrons in the Mott-Hubbard 
insulating phase. The chemical potential is /i = 1.08 and the 
occupation is n = Ud + rip — 3, with Ud = 1 + i^, rip — 2 — u 
and u — 0.14. Am denotes the Mott insulating gap. 



dynamical bath of the site, which is subject to the self- 
consistent condition. 



'-pd 
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(3) 



The p-electron Green function Gpp, is written in terms of 
the local self-energy Epp and the non interacting density 
of states as: 



Gppiiuin) 



de- 



Spp(i0„) 



(4) 



where e is the single particle energy, and Spp is obtained 
from the quantum impurity model self-energy S(ia;„) of 

Eq.mSS, 



Spp(z0„) — 



i0n + fi - ed - S(i0„) 



(5) 



To solve these equations we use two, in principle ex- 
act, numerical methods: Quantum Monte Carlo (QMC) 
[ist and Exact Diagonalization with the Density Matrix 
Renormalization Group (ED-DMRG) . The former is 
a finite temperature calculation and is exact in a statis- 
tical sense, while the latter is a T = method that relies 
on diagonalization of large finite clusters. The compari- 
son of the results obtained from these methods allows for 
a non-trivial benchmark of the numerical results. In our 
simulations we typically perform over 10^ Monte Carlo 
sweeps in order to minimize the statistical errors. 

In the main panel of Fig. [T] we show the DOS of the 
metallic state that results of lightly doping S = 3 — n 
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holes into the "parent" Mott insulator state (shown in 
the right inset). The main features of the metallic state 
are its mix- valence character and the presence of a quasi- 
particle peak at the Fermi en ergy that is flanked by lower 
and upper Hubbard bands [l6|. Within Fermi liquid 
theory, the key concept is that the one-particle excita- 
tions near the Fermi energy, the "quasiparticles" , are 
long lived entities. In consequence, the imaginary part 
of their self-energy should go to zero as ^ 0. More 
precisely, E(ict;„) « jw„a -I- b for small a;„, with a and 
b real constants. Quite surprisingly, however, the self- 
energy of the metallic peak of Fig. [1] does not have this 
property. In Fig. [2] we show the numerical results for 
various paramagnetic metallic states obtained by small 
doping of the Mott insulator. The main panel shows the 
large finite intercept of ImSpp that reveals the NFL na- 
ture of the states. Therefore, the metallic peak at the 
Fermi energy of Fig. [T] is not due to quasiparticles but 
rather to incoherent {i.e. short lived) excitations. In the 
lower inset of Fig. [2] we show the extrapolated value of 
ImEpp(tj — > 0) that is an estimate of the resistivity p{T). 
Contrary to what is expected for a metal, the resistivity 
increases with decreasing T. At higher doping, however, 
it crosses over to conventional metallic (ie, Fermi liquid) 
behavior. This feature is in qualitative agreement with 
experiments in various NFL heavy fermion compounds 
[l|. The coherence temperature scale Tcoh, below which 
the Fermi liquid behavior is obtained depends on 6. It 
is computed from the local spin susceptibility xioc aud 
is given by the x-axis intercept of l/xioc{T) T + Tcoh. 
Physically, it corresponds to the Kondo quenching of the 
local magnetic moments. We found that Tcoh decreases 
with decreasing doping, and becomes vanishingly small 
at low S. This feature is reminiscent of the phenomenon 
of exhaustion, where few conduction electrons have to 
screen a large number of spins. It has been extensively 
investigated in the Kondo lattice model [i^, [2l| and also 
in the PAM [H, [H, HI . The exhaustion is usually real- 
ized in the limit where the number of local electrons is 
restricted to unity, so to represent a spin at every lattice 
site, while the number of conduction electrons is taken 
vanishingly small. Here, in contrast, as Fig. [1] shows, we 
are in a more mix-valence situation where the number 
of p-carriers (holes) is not very small and never vanishes 
(it is ly — 0.14 at S — 0, and increases further with dop- 
ing). Nevertheless, in the present case also the exhaus- 
tion physics is in fact responsible for the vanishing Tcoh. 
One should realize, however, that the relevant number of 
" conduction electrons" is not given by the finite nominal 
p-hole occupation, but by S which is the number of doped 
carriers to the (mix-valance) Mott insulator state. 

Another important difference with respect to the stan- 
dard exhaustion scenario is the stability of the present 
NFL regime regarding magnetic ordering. The exhaus- 
tion regime is usually found to be unstable towards mag- 
netic phases i2^]. Here, in contrast we find that the NFL 




FIG. 2: Imaginary parts of the p— electron self-energy 
lmT,pp{iujn) in the NFL state, from QMC at T=l/128. The 
hole doping is 5 = 0.016 (squares), S = 0.049 (circles), 
S — 0.093 (triangles), and /i — 0.58,0.53 and 0.48 respec- 
tively. Solid line is the same quantity from ED-DMRG. The 
ED-DMRG data are obtained for large finite clusters of 40 
sites and are plotted down to the frequency of the lowest en- 
ergy pole. Inset: inverse scattering time ImEpp(tJ — > 0) as 
a function of temperature. Black symbols are for the same 
values of S as in the main panel. Open symbols are for higher 
doping S = 0.44 (circles) and 0.62 (squares) with /i — 0.23 
and 0.13 respectively. Note that the inverse scattering time 
goes to in the Fermi liquid state. 



regime is largely paramagnetic, except at very low dop- 
ings and temperatures, where the system becomes an an- 
tiferromagnetic metal (AFM). The reason for the mag- 
netic stability of the NFL regime can be argued to result 
from the magnetic frustration due to competing magnetic 
interactions. On one hand, neighboring localized elec- 
trons in the Mott state interact via the superexchange 
mechanism, that follows from the mapping of the narrow 
band to the Hubbard model. This interaction is antifer- 
romagnetic and is responsible for the AFM state that is 
found in the limit of very low hole doping. On the other 
hand, there is also a competing ferromagnetic interaction 
induced by hole doping. The doped holes need to delocal- 
ize to lower their kinetic energy, but they are subject to 
a strong on-site magnetic binding to the local moments 
[l^. Therefore, in order to hop, they need the local mo- 
ments of the neighbor sites to have the same magnetic 
orientation as the one in its current site. These com- 
peting magnetic interactions lead to unquenched local 
dynamical fluctuations of the magnetic moments down 
to the vanishingly small Tcoh, that provide the source for 
the NFL scattering. 

Fig. [3] condenses our results into a S-T phase diagram. 
The intensity plot shows the magnitude of the scatter- 
ing rate Iml]pp(ix' — > 0), which is a measure of the NFL 
character of the system. We observe that conventional 
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FIG. 3: Intensity plot of the scattering rate lm'Epp{uj 0) as 
a function of doping S and T. For visualization, the scattering 
rates are normalized to max{ImEpp(tj„)} at each {S, T). The 
line with square dots at low S is TnchI and gives the AFM 
phase boundary (The T=0 data-point is obtained from ED). 
The line with circle dots denotes the numerical estimate for 
the crossover scale Tcoh- It is obtained from the extrapolation 
of the low T behavior of l/xioc{T) (see text). 



metallic (ie, Fermi liquid) behavior occurs in two dark 
regions. The one at low 5 corresponds to the small AFM 
phase that we mentioned before, while the one at higher 5 
is the Fermi liquid state that is realized beneath Tcoh due 
to conventional Kondo screening. In between these two 
regimes we find the strong NFL behavior of the system. 

The normal metal behavior in the AFM phase is easy 
to understand. As soon as the magnetic moments become 
Neel ordered the doped carriers have no problem to form 
coherent electronic waves. Thus, a direct implication of 
this observation is that if one induces magnetic order, the 
NFL state should become normal. To test this hypothesis 
we have mapped out the phase diagram as function of the 
external magnetic field B and T at a small fixed doping. 
The results are shown in Fig. 3] which displays the inten- 
sity plot of the computed scattering rate. At S ^ we 
find a small dark normal metal region that corresponds to 
the AFM state discussed before. At large B the strong 
magnetic field induces ferromagnetic order of the local 
moments and, as expected, the metallic state is normal. 
Interestingly, in the intermediate B region, a NFL regime 
emerges which can be understood as resulting from the 
competition between the magnetic interactions. 

To conclude, we have studied the NFL regime that 
is realized upon doping a Mott insulator state in the 
Periodic Anderson model within the Dynamical Mean 
Field Theory. The NFL behavior is originated in the 
strong local magnetic scattering of the doped carriers 
by unquenched fiuctuating moments. This is in con- 
trast to other approaches where the NFL state results 
from scattering by spatial magnetic fluctuations. We find 
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FIG. 4: Intensity plot of the scattering rate ImEpp(aj 0) 
as a function of external magnetic field B and T, at fixed 
doping & = 0.01. For visualization, the scattering rates are 
normalized to max{ImEpp(a;„)} at each (B, T). 



that both, doping and magnetic field, allow to tune into 
the NFL regime, in qualitative agreement with some of 
the non-Fermi liquid phenomenology observed in heavy 
fermion systems. 
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